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Abstract. We investigate the differentiability of minimal average energy associated to the 
functionals Ssiu) = j^d jl^^P + eV(x,u) dx, using numerical and perturbative methods. 
We use the Sobolev gradient descent method as a numerical tool to compute solutions of 
the Euler-Lagrange equations with some periodicity conditions; this is the cell problem 
in homogenization. We use these solutions to determine the average minimal energy as a 
function of the slope. We also obtain a representation of the solutions to the Euler-Lagrange 
equations as a Lindstedt series in the perturbation parameter e, and use this to confirm our 
numerical results. Additionally, we prove convergence of the Lindstedt series. 



L Introduction 

, be periodic under integer translations. That 
is y(x -I- k,y -H = V(x,y) for all (k, /) e Z"' x Z, where (x,y) = (xi, . . .,Xd,y) e R'' x R. 
Furthermore, assume V is analytic. We consider the formal variational problem 

(1) Si.(u(x))^ I -\'^u(x)f- + eV(x,u(x))dx, 



[ ^IVm(x) 



where e is a small parameter, so that Se is a small perturbation of So(u) = J^j jIVm^ dx. We 
call u e Hl^iW') a minimizer (or minimal solution) of if for all </> e H^g^^iW^) 

r 1 , 1 , 

- 1 V(m + 0)1+ eV(x, M + 0) - - 1 VmP - eV(x, u) dx > 0. 

Usupp((p) ^ ^ 

Any minimizer of (1) must solve the Euler-Lagrange equation 

(2) - Am + eWyix, u) = 0, 

and by standard elliptic regularity theory will be at least as regular as V, so analytic in our 
case. 

Definition 1.1. Following [Mos86], we say a continuous function u is non-selfintersecting 
if the graph of u does not intersect integer translates of itself. That is, if m e C*'(R'', R) and 
V(fe, /) e Z"' X Z 

(3) u(x + k) + I - u(x) > 0, or < 0, or s 0, 
where the three alternatives are independent of x. 
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This is also referred to in the Hterature as the Birkhoff property. From a geometric view- 
point, this means that the graph of u projects into T''^' = R^'+'/Z''^' without intersecting 
itself, unless it coincides exactly. 

If M is non-selfintersecting, then there is a rotation vector, oj e W', associated to u such 
that 

(4) sup \u(x) - (L) ■ x\ < oo 



[Mos86]. A function ii satisfying (4) is called plane-like because its graph is at a bounded 
distance from the a hyperplane in R''^' with normal vector {a>, -1). 

We denote the set of all non-selfintersecting minimizers of Ss with rotation vector oj by 
MojiSe), and more briefly as Aioj, where dependence on is understood. As shown in 
[Mos86], Ma, is nonempty for all o) e W'. 

The rational dependency of w can play a role in the structure of Aia,- The case where 
o) = (w, - 1) is rationally independent was studied in [Ban87b], and the rationally dependent 
case in [Ban89]. 

We define the minimal average energy Ag : R'' — > R by 



where u e Al„, Br - {x e W' : \x\ < R). It was shown in [Sen91] that is wefl-defined, 
that is, the limit exists and is independent of the choice of minimizer u e Ai^i- Furthermore, 
Ag is convex, so that one-sided derivatives of A^ exist at each oj e R'', [Sen91]. In fact, for 
e = 0, any u e Aioj has the form u(x) = lj ■ x + a for some a e R (see [Mos86]). 
Thus, Aq(lij) = 1\cl>\^ is smooth. However, for typical V the differentiability of A,, breaks 
down when e > 0, and for large enough s the set of points where Ag is not differentiable 
will be dense in R'' (see [Ban87a]). Because Ag has one-sided derivatives the graph of 
Ag will have "corners" (i.e. the jump in the gradient of A^ in a direction e j) at points of 
nondiff'erentiability. In [Sen95], page 356 there is a formula for the one-sided directional 
derivative of A^ involving special types of minimizers of 5^, described in Section 3. 

In this paper, we present a numerical approach to computing solutions of (2) via a gra- 
dient descent method known as the Sobolev gradient, as explained in Section 2. We use 
this to compute A^ and Asia)) + D^e Aeiu), that is the jump in the gradient of A^ in the 
direction ej. We will sometimes refer to a jump in the gradient as a "corner". In Section 
3 we present a perturbation method for finding the minimizers of needed to apply the 
formula provided in [Sen95] to compute Ae(a;) + D^^Af^iio). 

Throughout the paper, we will deal with the case w e ^Z''. There are two good reasons 
for this. The first is that minimizers of (1) with irrational rotation vectors can be obtained 
as limits of sequences of minimizers with rational rotation vectors. The second reason 
is that for rational w, the partial diff'erential equations we need to solve are well-defined 
on the torus T''. That is, they have periodic boundary conditions. This allows the use of 
Fourier transforms in the numerical method. It also makes each equation in (16) of the 
form A0 - g, which can be solved for (p provided g has average zero. 

In this setting, for a fixed a> e ^Z"', the formal functional in (1) can be replaced by the 
reduced functional 



where z is A^Z''-periodic. Considering the limits of solutions un(x) = ojnx + Zn(x), cjn g 
j^Z'' as — » oo, we see that this problem is related to periodic homogenization. Rescaling 



(5) 




(6) 
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SO that (6) is defined on the unit cube, the process of letting — > oo is the same as letting 
the X dependence oscillate very rapidly. 



In this section we investigate numerically the size of the comers of Ag as s varies. Since 
our problem comes from a variational principle, steepest descent methods are a natural 
approach. We consider a fixed u e j^lf and seek solutions of (2) of the form u{x) - 
o) ■ X + z(x), where z is A^-periodic. To find such a function m, we solve 



for z, and set u(x) - lo-x+z(x). In this setting, we see that z is a critical point of the reduced 
variational problem 



We recall that the gradient of S^.n with respect to a Hilbert space, H, is the unique element 
ofh€H such that DSe,N{z)T] — {g, rj)H for all rj e H. For instance, the L^-gradient of Se,N at 
z is the unique element of L? = which we will write as VQSs,Niz), such that DSe,N{u)ri - 
(Vo5£,/v(z), T])i2. Integrating by parts, we have, DSe,N(z)r] = (-Az + eVyix, tj- x + z), rfJii. 

Considering gradients with respect to other inner products has been a fruitful endeavor 
(see [NeulO]). In our particular case, by considering the gradient of Se.N in for (3 e (0, 1] 
we avoid the stiffness of the problem that appears in the case. This is explained in more 
detail in the remark following the derivation of V^5£_Af(z). 

The standard inner product on H^{NT') is {u, v}fffi = ((/ - A'fu, v)//o. For y > we 
define the inner product (m, v)^ = {(yl - A'fu, v)//o, which determines a norm on that is 
equivalent to the standard norm. For more details on this see the introduction of [BLVll]. 
In that paper, the authors show that the descent equation d,u - -VpSg^Niu) satisfies a 
comparison principle, and preserves the class of non-selfintersecting functions. In this 
way, they obtain critical points of Se,N with rotation vector cj by choosing initial condition 
u{x, Q) = (L) ■ X, equivalently, z{x, 0) = u(x, 0) - lj ■ x = Q. 

We note that different choices of y result in different inner products, and therefore differ- 
ent gradients of Si;^^. Thus, '^pSg^u) depends not only on the choice of Hilbert space, but 
also on the choice of inner product on that space. We calculate the //'^-gradient as follows: 



2. Numerics 



(7) 



Az = eVyix, oj-x + z), z{x + Nk) = z(x), Vfc € if. 




where NT' - W'/NZ''. The Frechet derivative of Se,N at z, applied to tj e is 





{(yl - AfiyI - A)-^(-Az + sVy{x, u-x + z)), 77)^^ 
((y/ - A)"^(yz - Az - yz + sVy{x, <jj- x + z)), z?)^ 
((y/ - A)'-^ z - (y/ - A)-^(yz - sVyix, u-x + z)), /?) . 



Thus, our steepest descent equation in H^^, d,z - -^/sSe^z), becomes 
(8) d,z = - (y/ - A) z + (y/ - A)"^ (yz - sVyix, oj-x + z)). 



4 



TIMOTHY BLASS AND RAFAEL DE LA LLAVE 



z '■ [0, A^]'' R with periodic boundary conditions. 

If the solution t) of (8) approaches a critical point of Se,N, that is zix, t) — > zdx) as 
f — > oo, then Zc will solve (y + A)'"'' Zc - {j + ^^^{jZc - Vy(x, w ■ x + Zc)), which reduces to 
-AZf + eV'y(x, w ■ jc + Zc) = 0. Then = w ■ x + z^, will solve (2) and sup^ \u{x) -u- x\ < oo. 

Remark 1. The stiffness of (8) is significantly reduced for values of yS ^ 1. For in- 
stance, the gradient descent equation would be 5,m = Am - eVyix, u). If we set G{k) - 
T{-eVy(x, u)}(k), then the descent equation in the Fourier domain becomes the system 
dfUii = \k[^Uk + G{k) for each k. As \k\ grows, these equations become increasingly stiff (the 
eigenvalues of the linearization have a large spread). However, the H^^ gradient descent 
equation in the Fourier domain becomes dtUk - {y + |A;p)'"^Mi + (y + \k'^Y^G{k). So for 
yS ~ 1 the stiffness is greatly reduced. 

2. 1 . Implementation of Sobolev gradient descent. We fix w e ^Z'' and seek z : [0, oo) x 
[0,A^]'' — » R solving (7) with periodic boundary conditions. Thus, (7) can be rewritten 
in the frequency domain via the Fourier transform. The main benefit of this is that the 
pseudo-differential operators in (7) simplify greatly in the frequency domain. However, the 
composition Vy{x, o)- x + z{x)) is complicated by the Fourier transform. We take advantage 
of the simplicity of the operators in the frequency domain in our numerical scheme, and 
pay the price of computing an inverse fast Fourier transform at each time-step. This is 
explained in the remark following equation (11). 

We now develop the details of the implementation for d - 2. We know that z{x, t) - 
u{x, t) - w ■ X will be A^-periodic if we have Nu e 7? and if z{x, 0) = uix, 0) - w ■ x is an 
A^-periodic function. So we set zq(x) - z{x, 0) — 0, and consider the equation satisfied by 
z, the Fourier transform of z- We get 

(9) z,{t,^) = -(r + \^\^)'-Pz{t,o + (r + |^P)-^(rz(f,^) - nVy{x,co- x + z)](t,^)), 

where we have also used Ti-] to denote the Fourier transform. We now choose a time 
step Af and set t„ - nts.t. We also break up the domain [0, A^]^ into discrete points 
and represent z(f„,x) as an m x m array ZniUj)- So, representing the discrete fast Fourier 
transform as fft[-] (and with z„ - fft[z„]) equation (3) is approximated as 

(10) ^jiii^ ^ + |^h'-^i„+i + (r + I^h"^(r2«+i - frt[y,,(x, « ■ x + z„)]). 

This is a quasi-implicit method, since nearly all of the right-hand-side is evaluated at the 
later time t„+\. Only the nonlinear term is evaluated at time f„, so we can easily solve for 

Zn+i- 

. . ^ Ui, j) - Atjy + + ^2{jfrm[Vyix, oj ■ xji, j) + z„(/, ;))] 

i+Af(r + ^i(02 + ft(;)2)i-^-rAf(r + ^i(02 + ^2(;W ' 

Remark 2. After the «-th step, we have z„. In order to compute Vy{x, oj-x+Zn) we apply the 
inverse fast Fourier transform to get z„ = ifft[z„], which requires order log(OT) operations 
because z„ is an m x m array (or a vector of length m^). With z„ computed, we evaluate 
Vy(x, o) ■ x + Zn), requiring operations, and then we transform it with FFT, again costing 
0{m^ log(m)) operations. 

To compute z„+i, we now only need to perform component-wise multiplication of the 
arrays in (11), requiring 0(m^) operations. That is, we multiply the (/, j)-component of 
fft[y,.(x, w ■ x(i,j) + Zn{i,j))] by the (/, y)-component of At(y + ^i(i)^ + ^2(j)^T^ for each 
1 < i,j < m. We then add to it the (/, y)-component of z„ii,j) and then divide by the 
(/, j)-component of 1 + At(y + ^,(/)2 + ^zUf y-" - jMy + ^i(0' + ^jfT"- 

If z„ and fft[yj,(x, a> ■ x + z„)] were represented as vectors of length m^, then this same 
procedure would amount to multiplying fft[y,,(x, a> ■ x + z,,)] by the x diagonal matrix 
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representing Af(y + ^[ + ^f) ^, adding z„, and then multiplying the result by the x 
diagonal matrix representing (1 + At(j + + ^2)' ~ yMj + + ^2)"^) ' ■ 

Remark 3. This method would also work in a setting where the gradient part of the en- 
ergy functional S'g is replaced by the fractional laplacian. We could use all of the same 
techniques above for solving the gradient descent for 

(12) 5f(z)= r z(-A)^z + sV(x,aj-x + z)dx, 5 e (0, 1). 

The critical points of (12) will solve the Euler-Lagrange equation 

-(-A)*z = sVy(x,cL>-x + z). 

Using the metric on given by (v, w) = {{y + (-A)*y^v', w}[2, we arrive at the descent 
equation 

(13) d,z = -(r + (-A)^)'-^z + (r + i-Afy/^iyz - Vyix, OJ-X + z)). 

Because z is periodic and the operator (-A)* is diagonal in the Fourier coefficients, 
implementing (13) numerically is the same as described above, except with powers of 
in place of |^|^. 

3. Perturbation method for computing minimizers and the jumps in DAe{oj) 

3.1. Foliations and laminations of minimizers. Following [Sen95], we define = 
{(k, I) e Z'' X Z : oj ■ (k, I) - 0), where oj - (oj,-!), and M{a)) c Moj, the set of maximally 
periodic u e by 

M(ci)) - {u e Moj : u(x + k) + I - u(x), V(^, e F^) . 

For each a> e the set A1(a<) is closed and totally ordered. The closedness is a classical 
argument (see [Mor73]). The total order, meaning that if u, v e Ai{a>) then either m > v or 
M < V or M s V, is a consequence of the maximum principle for elliptic partial diff'erential 
equations. Let x - (xi, . . . , xj), and we write (x, Xd+i) G Ef''^^. Let oj G W^, then the total 
order of Ai{(l)) means that for a given (x, x^+i) e R''^' there is at most one u e Ai((l)) such 
that Xd+i = u(x). That is, each point in R''^' belongs to the graph of at most one u e Ai(a)). 
If for all (x, x^+i) e R"'^' there exists u e M{(jj) with Xd+\ = u(x), then we say A1(a)) is a 
foliation of R''^'. Because of the non-selfintersection property (3), such a foliation projects 
into T''+'. It can happen that there are points (x, x^+i) e R''^' for which there does not exist 
any u e Ai{d)) with x^+i - u{x). In this case we say A1(w) is a lamination of R''"^^ (and 
projects to a lamination of T''^^). For this reason, a lamination is sometimes referred to as 
a "foliation with gaps". 

If a> is rationally dependent, and if Ai(a>) defines a lamination, then there are minimizers 
u € Moj whose graphs lie in the gaps of Aiicu)), [Ban89]. In addition, if we choose a 
direction 

yS G spanjilfe eZ'' ■.Lo-keZ}n 5''"' 

then there is a m e Aia, such that u is asymptotic to some e A1(w) in the direction /3 and 
asymptotic to some e Ai{(h) in the direction -/3. For details of the asymptotic behavior, 
see [Sen95], page 350. Such a m is said to be heterocUnic in the direction fi. A formula 
for the one-sided directional derivative of Ag at a point w is given on page 356 of [Sen95], 
and involves integrating the action over the gaps defined by the elements of Micj) and the 
heteroclinics between them. We will use perturbation methods to calculate the gap borders 
and the heteroclinics lying in the gaps. 
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3.2. Lindstedt series for solutions. We seek plane-like solutions Ugix) of 
(14) Au-eVy(x,u) 

that can be expanded as Ueix) - UQ + eux+ s^U2 + The series 2j>o s-'uj shall be referred 

to as the Lindstedt series for the solution m^. Substituting the series into equation (14) and 
matching powers we arrive at the following equations for each order of s 





Am() 




Am] 


a": 


Am2 


(15) : 


Am3 




Am4 


In order that 


be 



1 2 

-VyyyiX, Uq)U ^(x) 
VyyyiX, Uq)Ui{x)U2{x) + 



each i > I- Using the notation to refer to the j-th coefficient of the power series in e, 
we write 

Vy{X, Ue) - Vy(X, UO+SUI+ ...) - '^[Vy(X, Mj;)] jB''. 

We will also write u^-i for the first j terms in the Lindstedt series: m"^-' = mq + ■ ■ • + s-'^^uj^i. 
The j-th order equation in the list (15) has the form 

(16) ej : Auj = [Vy(x, m)],-i = [Vy(x, u<j)] j^,. 

The zeroth order equation is satisfied by any affine function, so we take uq - a> ■ x -¥ a, 
and at this point we are free to choose a as we like. To solve the j-th order equation, 
we must have that jj.j[Vy(x,u'^-')]j-idx - 0. This compatibility condition is what forces 
specific choices of a once to has been fixed. 

We also note that the solution of (16) is determined only up to an additive constant, 
which will be chosen so that the equation of the following step has a solution. That is, the 
average of uj is chosen so that equation for uj+i is solvable. 



3.3. Existence of the Lindstedt series to all orders. We consider co e j^Z'' fixed, and 
seek Ug solving (14) such that u^ix) - a> ■ x e L^iNT^). We must have mo(x) = w ■ jc + a to 
solve Amo - 0. The choice of a is free, and each value of a e [0, 1) will result in a different 
set of equations for the uj with j > I. In this section, we show that there are at least two 
choices of a e [0, 1) such that (16) has a solution for each /' > 0. 

Lemma 3.1. For fixed to e ^Z"', there are at least two choices of a e [0, 1) such that 
j^T<i Vyix, oj- x + a)dx^Q. 



COMPUTING THE MINIMAL AVERAGE ENERGY 



7 



Proof. Define Oi : R — > R by ^\{a) - Vy(x, to - x + a) dx. Then Oi is continuous, and 
we have 



I ^i(a) da = I I Vy(x, a) ■ x + a) dxda 
Jo Jo Jnti 



Jtrc Jo 



d 

— V(x, tu ■ X + a) dadx 
da 



= I V(x, a> ■ x) — V(x, a> ■ X + l)dx = 0. 

Thus, must have a zero in [0, 1), and since <l>i(Qf + 1) = 'i>i{a), by the periodicity of V, 
we know that Oi must have at least two zeros. □ 

For any such choice of a, Vy(x, mq) dx - 0, and therefore there exists a family of 
periodic solutions, ui(x) - m*(jc) + A, of Amj - y,,(x, mq), differing only by an additive 
constant. We will write u* for the member of the family with average zero. 

Theorem 3.2. Let uq - oj ■ x + a, with u e j^Z^ and (b\{a) -Q. If 

(17) r yj.,(x,Mo)t/x^O 

Jnt' 

then each equation Auj — ^Vy{x, M*-')j . ^ has a solution for all j > 1. 

Proof. From Lemma 3.1, we have a family of solutions mi = m* + A, of Ami - Vy(x, uq). If 
we set 

then, [yy(x, M'^^)j^ c/x = J^^^ y,,y(x, mo)(Mi + = 0, so the equation Au2 = 

yvy(x, mo)mi is solvable for U2 - U2 + A when we choose U]{x) - u*^(x) + A'-^K To con- 
tinue inductively, we note that for each j > I, 

[VyiX, M^Oj^. j - Vyy(X, Mo)Mj-1 + Rj{Ul,. . . , U j^l) ■ 

Suppose that we have solutions of (16) for j - 1, . . . , n, and the constants A^^\ . . . , /i*""^^ 
have been selected so that Uj(x) - u*(x) + A^~'\ We choose 

^(«) _ ivT" ^yy'^^' +Rn{UU---,Un-l)dx 

f^,Vyy(X,U0)dX 

so that 

I \Vy(x,u^"^^)] dx^ f Vyy(x,uo)(ul+A''"'>) + R„(uu...,u„^i)dx^O. 
Jnti " Jw 

Thus Am„+i = [y^.(jc,M<«+i)]^^ has a family of solutions, u„+i(x) - m*^j(x) + A, completing 
the induction. □ 

3.4. Convergence of the Lindstedt series. We use a Newton method to produce a se- 
quence of functions f/„(x, s) that are analytic in e, and converge uniformly to a solution of 
-Aue + eyv(jt;, Me) - for e e ^^(0) c C for small enough 6 > 0. Thus, we produce an 
e-analytic function, Ug, solving (14) for e e Bs(Q), so the Taylor series of Ug must coin- 
cide with the Lindstedt series, proving the convergence of the Lindstedt series. For similar 
convergence results, but for the case of Diophantine frequencies, see [CF96]. 
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Lemma 3.3. Let uj 6 i/"'^^(T'') and e e C and define u{s,x) = Yij£N^'''^j(^)' where the 
series is convergent in H""^^ for \e\ < r, for some r > and m > d/2. Define F : C — » 
i/'"(T'') by 

F(e; u) — Au(e, x) + eV(x, u(e, x)) 
Then the derivative of F with respect to e is 

DeF(s; u) — ADsu(e, x) + V(x, u(e, x)) + sVyix, u(e, x))Deu(e, x). 
Proof. We define 

G(e; m) — ADi;U(e, x) + V{x, u(e, x)) + EVy(x, u(e, x))Di.u{e, x) 
H{e; u) — ADIu(e, x) + 2Vy(x, u(e, x))Deu(e, x) 



We have 



+ BVy(x, u{e, x))D\u{e, x) + sVyyix, u{e , X)) (DeUiE, x)}^ . 



I H{e + cnh; m)t/i dcr — G(e + rh; u) — G{b; u) 
Jo 

JG(e + rh; u)h dr = F(e + h; u) — F(e; u) 




so that 



and 



; u)h =1 I 

Jo Jo 



F(e + h;u) - F(e; u) — G(e; u)h — I I H{e + cnh; ujrh da dr 



\\H(e;u)\\h„. < \\Dlu(E,x)\\^,„^,+2\Vy\^, \\DMs,x)\\h,.. 

+ \e\ |y4„ \\dIu{e,x)\\jj,„ + \e\ |y,4„ ||(Z)^m(£,x))2||^„, . 

From the Gagliardo-Nirenberg inequality [Nir59], we have that (DgU{E, xyf e //'" because 
m > d/2, and that there is a constant depending on m and d, such that ||(D£M(e, x))^||^„, < 
C(m, d) \\Deu{e, x)\^fj,„. Hence, 

\\F{E + h;u)-F{E\u)-G{B\u)h\\H„ < f f \\H(e + a-rh;u)\\fi„T\h\^ do-dT 

Jo Jo 

<C(m,d,\\u\U2AV\c2,£)\hf. 
Thus, F is difFerentiable and D^FiE; u) - G(b; u) □ 

We define /^^(m) = -Am + sVyix, u). It was shown in the Section 3.3 that for any fixed 
M e N, we can solve the first M equations from (15) for m^*^, and that m^'^ will solve (14) 
up to order e'^. That is, F^iu";^) = 0{e'^). We set t/o = for a sufficiently large M to 
be determined later, and for 7 > 1 we define 

(18) f/„+i = f/„ - DF,{U„)-'F{U„). 

Let m > d/2 he fixed, and note that for any M > 0, u^'^ e by the regularity theory for 
elliptic PDEs. If t/„(x, e) is analytic in e and is jjj jjjgjj Fe{U„) - -AU„+EVy(x, U„) 
is analytic in e and //'" in x by the result in Lemma 3.3. We have OFsiUn)^ - -A77 + 
sVyyix, Un)ri, and we need to consider carefully the behavior of DFe(U„y^ . To simplify 
notation, we define L'l : H"'+^(NT') H"'{NT'') as 

L: = DF,(U„) = -A + EVyyiX, U„). 

L'l is a small perturbation of -A : _^ jsJqw, -A maps the codimension one 

subspace H"'^^/R of its domain to the codimension one subspace //"/R of its range in a 
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bounded, invertible way. But it has the simple eigenvalue A - 0, with eigenspace spanned 
by the constant functions. 

Lemma 3.4. Let Pq ; H'" — > H'" /M. be the orthogonal projection onto H"^ functions with 
zero average. That is, iff e H"\ then Pq/ = / - fi„, fdx. Let -Aq : H'"^^/R H'"/R be 
the restriction of the laplacian, and let Y be the image of H'"'''^ /R. under — Aq + eVyyix, U„). 
Suppose there are constants q,ci,6Q > such that 

(19) (sVyyix, Un),g)H'« > cilel^ e y\ = 1 
for \e\ < (Jo. Then ||(L")"'||x(ff") < ce^'', for some c> 0. 

Proof Let g = -Ao + PosVyy{x, U„), then Q : H"'+^/R //'"/R, and for small s will have 
a bounded inverse, g^' : H"'/R H'"^^/R. 

We have that -Ao + sVyyix, U„) = -Ao + PQsVyy{x, U„) + P^sVyy(x, U„) maps H"'^^ jR 
into i/"' and is a small perturbation of Q, and there is a C2 > such that || - AqT] + 
eVyy(x, Un)ii\\H"' ^ C2||77||h"'+2 for all rj e H"''^^/R. Thus, F is a codimension one linear space 
isomorphic to 7/'"/R lying inside H'", and - {g e H'" : j^^^j(-AoT]+sVyy(x, U„)rj)g dx - 
0,V/7 e H'"^^ IR] is one dimensional. 

The condition (19) implies that the image of constant functions under L" does not lie 
entirely in Y, but has a component in the Y^ direction. Thus the image 

under L" IS 

H"\ and for small e, L" will be invertible. 

Let Py and Py denote the orthogonal projections of H'" onto Y and Y^. We let g e Y-'- 
be a unit vector such that Py^ - {^,g)H'"g for ^ e //"', which is possible because Y^ is 
one dimensional. Let ^ e H"', so ^ = L"ri for some rj e H'"'^^. We write rj - rji + rjo with 
T]i e H"'+^/R, and t]o & R, and ^ ^^y + ^i. with ^y e Y and e F^. Writing ^ in terms of 
771, 7/0 we have 

At; + sVyyix, Un)r] = Aqt/i + sVyy(x, U„)rii + sVyyix, Un)r]o 

= A077I + sVyyix, U„)ril + PySVyyiX, U„)T]Q + PySVyyiX, f/n)/?!). 

The term A^rji + sVyy{x, U„)ri\ + PYsVyyix, U„)riQ e Y, so the component of L"ri in F-"" is 
^± = {L"r],g)H'''g = <eyvv(jc, t/«)?7o,^)ff"g- Thus, 770 = ^±«eVvv(-«, t/„),^)/f».)"' and t/j is 
given by 

,,.{A..sVy.ix,U„)) 
Hence, ||(L«)-'^||ff„.2 < -^||^J|^,„ + -L (||^y||^„, + e||y„,||H." ^ll^JIff-) < c|er?||^||„».. n 

Recall that (i") ' is a compact operator by the regularity theory for elliptic PDE. In 
particular, the eigenvalues of L" are isolated from the spectrum, and if A„ is an eigenvalue 
of L", then the resolvent of L", written as R{L", = (L" - (ly^, has the representation 

(20) r(l:, = V(^ - O'Qi"' + 1 — -Pn, 

where Q„ is bounded, and P„ is the spectral projection on the the A„ eigenspace: 

P„^-^.j'^R((,L"^)d(, 

and F is a closed curve enclosing A„ but no other point of the spectrum (see [Kat66]). The 
principal eigenvalue, Aois), of L" is simple because L" is an elliptic operator, [Eva98]. This 
means that Ao{s) is analytic in a neighborhood of e = 0, [Kat66]. In the iteration process. 
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{LIY^ = R{LI,0 = 27=0 'InQn^' + f^B' will act on F^{U„). At each step, we need the 
function U„ to be analytic in e, so we need {L"y^ Fi;(U„) to be analytic. 

Proposition 1. Assuming condition (19) holds, there is a choice of M e M such that, 
if Uq — u^^, then U„ will be analytic in s for all n > 0, and the Newton method (18) 
converges uniformly in s in a neighborhood of s — 0. 

Proof In the iteration process, {L'l)-^ = R{L1, = S™,, A^QC^ + j-P„, will act on F,iU„). 
As explained in the previous paragraph, Ao(s) is analytic in a neighborhood of e = 0, and 
will have a zero of order p e N at e = 0. From the result of Lemma 3.4 we may assume 
p<q. 

We will use induction on n to show the analyticity of U„+]. Fsiu^'^) has a zero of order 
M at £ = by construction of u^'^, and we take M > 2q. By the expansion in (20) and the 
result from Lemma 3.3, for ^ = and A„ = Ao(s), (L^l)^^ F ^iu"^^) has a zero of order M - p 
if M > /? or a pole of order p - M if p > M. Thus, taking M > 2^ is more than enough 
to ensure U\ - m^*^ + (Lg)"'Fe(M^'*^) is analytic in e, since q > m. From (18), FgiUi) = 
D^F^{Uq){DF^{Uq)-^F^{Uq))^ + . . ., and since D^FJUo) = sVyyy(x, Uq) is bounded, we 
have F,(f/i)||ff» < CN(||(L°)-i||x(ff")||F,(f/o)||H™)'. Hence, \\f's{Ui)\\h«^ < Qsl'-^"^^^ 
by Lemma 3.4 and because ||/^£(t/o)llH"' ^ kl'^- We chose M such that M - q > q, so 
\\Fe(Ui)\\H'« < C\s\^' where Mi = 1 + 2(M - q) > 2q. We also have \\Ui - t/ollff" < c\s\i, 
so 

Ll^-A + sVyyiUi) = -A + sVyyiUo + f/i " Uq) = + Ois") . 

Hence, L' is a small perturbation of L^, and Ai{e) will have a zero of order p. Thus, we 
have established the first step of the induction on n, because Ui is analytic in s, Fs(Ui) has 
a zero of order Mi > 2q, and the principal eigenvalue has a zero of order p < q at s - 0. 

Now assume that U„ is analytic in s and Fg{U„) has a zero of order M„ > 2g' at e = 
0. FJUn^i) = D^F,(Un)(DF,(U„)-'F,(U„))^ + ... and \\D^FJU„)\\nH'n < e||VvvvllL» is 
independent of «. So 

(21) ||F,(f/„+i)||ff„ < Ce||(DF,(f/„)-'F,(f/„))||^,„ < C£(e*^"-'')2. 

From (21) we have that U„+i is analytic, and FsiUn+i) has a zero of order at least 2q at 
e = 0. Just as for the n - Q case, \\U„+i - U„\\h'» < c\sf so An+i(e) has a zero of order p < q 
at e = 0, and the induction is complete. 

To prove uniform convergence, we have from (21) the recurrence formula: 

\\Fs(U„)\y« < (e-''\\FJU„^i)\M\ 

provided we choose e small enough (independently of n) so that Cs < 1. Therefore, 

\\Fe(U„)\\H»' < \e\-^''^^"-^^\\Fe(UQ)f^,„ < c\sf(M-2g)+2q_ n^j^^e, with M > 2q, the error on 
FsiUn) is bounded by cs^" , c independent of n. □ 

3.5. Connecting orbits and corners of the energy. The connecting orbits, or heterochnic 
orbits, exist only when Micb) fails to produce a foliation of T''^', as described at the be- 
ginning of Section 3. Each u e Ai{d)) satisfies Am = sVy(x, u) because it is a minimizer of 
S E. The order- zero approximation to u has the form uq - oj ■ x + a, for a e [0, 1). This is a 
continuous family, whose graphs foliate T''+'. 

From Lemma 3.1 we know that if V satisfies the twist condition (17) then there are at 
least two choices of a e [0, 1) such that the Lindstedt equations (15) can be solved to all 
orders. Recall that the function Oi, defined in Lemma 3.1, has at least two zeros in [0, 1). 
If the zero set of Oi is not all of [0, 1), then a choice of a must be made in order to solve 
Ami = Vyix^, cl>- x + a) for periodic mi. Thus, the order-e approximations, given by mq + sui 
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are not a continuous family indexed by a e [0, 1), but rather by a strict subset of [0, 1). The 
graphs of the functions in t his family no longer foliate T"'^' . 

If Oi(Qr) = 0, no gaps appear in the approximation up to first order in s. In this case, 
we can find periodic mi solving Amj - y,.(x, w ■ x + a) for arbitrary a, and the order-e 
approximations do foliate T'^'. We can then move on to try to solve Au2 = [Vy{x,u'^^)]i 
for periodic U2. With a free, we define : [0, 1) — > R by 02(0;) - J^[Vy(x,u'^^)]idx. 
Now either <I>2(ci') = or there is some a € [0, 1) with <1>2(q') 9^ 0. In the latter case, we 
must make a choice of a such that ^2(0:) - 0. In the former case, we move on to find M3. 
This can continue as long as ^j(a) - J^^j[Vyix,u'^-^)]j^idx is identically zero. Otherwise, 
a choice of a must be fixed, and the foliation breaks down. 

If at the j-th step we find that Ojia) + for some value of a then we make the ansatz 

that the heteroclinic orbit will be of the form Uh(x) - u) ■ x + a(e-'^^x) + eu\(x) + 

And a(e'^^x) ^ a-t as x — > +00 where a± satisfy <bj{ai^ - 0. Then Aa = 0{e') and 
the J-th order equation will be Aa + Am^ - [yy(x, M^')]j_i. Thus, a will solve a PDE 
of the form Aa = f{x, a) with boundary conditions a{e-'^^x) — > a-t as x — > +00, where 
J^jf(x,a)dx + 0. That is, / is the term that keeps [Vj^jc, M*^-')]j_i from having a zero 
average for any a. This is carried out in detail in Section 4 for a specific example. 

3.6. Energies involving the fractional Laplacian. The existence of minimizers in the 
case (jj € ^Z'' has been established for energy functionals involving fractional powers of 
the laplacian, [LV09], [BLVll]. As described in the remark at the end of Section 2.1, the 
Euler-Lagrange equation for the functional 

Sl{u)- I u{-t^^u + eV{x,u)dx, 5 e (0, 1), u - u ■ x + z{x), 

is 

—{—Sfu - Yy{x, u), u - (X) ■ X + z{_x). 

Much of what has been described so far regarding solutions of Am - Vy{x, u) carries over 
to this case. However, the analogous properties of the associated minimal average energy 
A* that are presented for <5 = 1 in [Sen91] and [Sen95] need to be established if one desires 
to investigate the differentiabiUty properties of in this case of non-local energy. This 
remains an interesting challenge. 



4. An example of the Lindstedt series 

Consider the potential V : R"'^' — > R given by V{x,y) - sin(27rA: ■ x)cos(27ry), with 
k e if. For a fixed w e R'', We compute the jumps in the derivative of the average energy 
functional, Af;(u)) defined in (5), using asymptotic expansions of the connecting orbits. 

Let Ue{x) solve Am^ - eVix, Ue). Writing - YjJ we see that Amq - 0, and thus 
uo = to ■ X + a, a> e W^, a e R. We are considering w fixed, but the choice of a is free and 
we have a family of solutions parametrized by a e R. 

However, in order to calculate ui we must solve Ami = Vy(x, mq) for a periodic function 
Ml. The average of Vy{x, mq) depends on a>, so in this example we will choose w = ^ so 
that the average of Vy{x, mq) is nonzero for some choices of a. Thus, the requirement that 
Vy{x,uo) = -27r sin(27rA: ■ x)sin(27rA: ■ x + 2na) has zero average forces cos(27ra') = 0. We 
say that Ue has a first order resonance if a is restricted in solving for the order e' term in 
the Lindstedt series. 
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We can calculate the first two terms in the Lindstedt series for each admissible value of 
a - 1/4, 3/4. We have Ami = +7t siniAnkx) so that the two solutions are: 

ul^kx+l - / siniAnk ■ x) + 0{s^) 
4 IottIkI^ 



1 e 



u"! - kx + - + r- sin(47r^ ■ x) + O(s^). 

4 167r|A:|'^ 

The superscripts c and m indicate that the solution for cc = 1/4 is a minimizer of the energy 
functional 

1 r 1 , 

Ae(u) - lim — 7 I -\Vu\ + sV(x,u)dx, 

N^oc N'' J[o,iv]<' 2 

and the solution for a = 3/4 is a critical point, but not a minimizer We restrict our attention 
to minimizers, where Ag can be considered as a function of the rotation vector w. 

To compute the jump in the gradient of at w = ^ using the formula in [Sen95] we 
need to find minimizers of Ag that are heteroclinic between u"' and m™ + 1 . We search for an 
asymptotic solution to Am = sVy(x, u) that has the form u'lix) - kx + a{ -sfex) + eui(x) + o(s). 

The order s' equation is 

Aa + Ami = -2?: sin(27rA:x) sinilnkx + 2jia) 

(22) 

- -n cos(27ra') + n cos(2;rQ') cos(4;7rA:x) - n sin(27rQ') sin(4;7rA:x). 

We want to choose a so that it is essentially one-dimensional (i.e. for some t] e R"', we want 
a to be a function of 77 ■ x). It should also satisfy Aa - -n cosilna) and we can eliminate 
those two terms from the equation above. At this point, any choice of rj is reasonable, and 
in the end we will chose rj to be in the direction in which we want to differentiate Asicj) 
(thus, 77 will typically be a standard basis vector). 
Letting fj denote ^^rj, we have 

1 r 3 

a(z) - - arctan(sinh( yljTfj ■ z)) H — 

7T 4 

Recall that in the expression for m, a is evaluated at z = ■\/bx. Expanding sin(2;Ta( s/e' x)) 
and cos(2:7ra( ■\/s' x)) in Taylor series, for this choice of a, we have 

sin(2.a( VFx)) = - 1 . 2( V2^./) • xf - 1)( ^nf^ ■ x)U . . . 
= -1 + 0(e) 

cos(27ra( ^fF x)) - 2( ^^lenf] ■ x) - ^( ^^Isnf] ■ xf' + . . . 

= o(V^). 

Then equation (22) becomes Amj = tt sin(47r^ji:) + 0{ VF ) and we get for mJI the expression: 

M*(x) - k - X + a( Ve~ x) + eui(x) + 0{s^^^) 

1 I — 3 e 

- k ■ X -\ — arctan(sinh( y2s ni] ■ x)) -\ r- sin(47rA:x) + o(e). 

n 4 16;r|^|'^ 

To be precise, we should write m^ ^ to show the dependence of m^ on 77. It is important 
to notice that m^ ^ is heteroclinic from m™ + 1 to m"' in the direction -rj. Similarly, u'l _^ is 
heteroclinic from m™ + 1 to m"' in the direction tj, and from m™ to m™ + 1 in the direction -77. 
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4. 1 . Computing the gradient of Asioj). Here we restrict 77 to be a standard basis vector 
e j. To ease notation we will write M(x) - m™(x) and He.{x) - u'l^ix). 

A formula for the derivative A^Cw) is found in [Sen95], and in our case becomes 



r r°° 1 1 

D.A^iu) = -|VM|2 + eV{x, M) - -|V//,/ - eV{x, H,) dxj dx''-'' 

' J[0,1]''-' J-oa 2 2 ' ' 

r r°° 1 1 

D_,,A,(w)= -\VH_,.f +sV(x,H_,^)- -\VMf -EV(x,M)dxjdx' 

Jm.ii''-' J-00 ^ ^ 



''[0,1]' 

Thus, the derivative of A£(w) in the direction Cj is difference in the energies of the minimizer 
defining the top of each connected component of gap and the minimizing heteroclinic in 
the direction Cj. Note that |V(M + 1)P + sYix, M + 1) = |VMp + ey(x, M). 

4.2. Two dimensional case. We now focus on the two dimensional case so that we can 
compare with the numerical computations, which were done for the two dimensional case. 
The computations in higher dimensions are impractical for us at the moment. 

In our example, V//,, - k + rj yflej cosh( y/lenrjx) - ke cos(4nkx) / 4-\k\^ . If we select 
77 = ei and compute Dg^Sik) we have 

1 T 1 o k\ \1e e e m 
-|Vi/„|2 ^ |^|2 + + ^ cos(47r^x) + 0{E^I^). 

2 2 cosh( V^eTrjci) cosh^( V2b ttjci) 4 

The heteroclinic solution has the same expression as above, with a sign change 
on the second term. The \k\^ terms from |V//±(.J and |VM| = + (9(e) will cancel 
when computing D^i,^Ae{u). When computing the jump in the derivative, that is the sum 
De,As{a)) + Ae(c<j), the terms ± — — ^{\\ cancel, and we are left with two terms 

cosh( V2e nx\ ) 

of the form — , " — — . These two terms then contribute for a total of 

cosh"( yle tixi) 

dx\ 



cosh^( V2e~ TTXi) 



If we integrate over all of R then this integral is 



%J — cx 



Is 2V2 ,„ 
dx\ - e ' . 



cosh^( V2e' nx\) 



Computing the contributions from the potentials, eV{x, M) - eV(x, Hj-i,^)dxj dx'' 

is not easy analytically, but numerically we find that they also yield ^^g'^^. Thus, the 
jump in the derivative is Di,,Ae(k) + D^e^Aeik) - i^g^/^ ^ l.ge'^^, which agrees well with 
the numerical computations. 

4.3. Comparison with results from the numerical computations. We use the Sobolev 
gradient method to compute the minimizers for several values of e, where each has 
rotation vector co. We then repeat this computation for the same values of e, but for rotation 
vectors oj + Acjej. With these minimizers we can compute As{oj) for each value of s 
and each rotation vector Taking differences over the w-variable, we can approximate the 
derivative of Ae(w) with respect to lj. 

D. AJoj) K — ^ . 

The plots in Figures 1 and 2 are for three cases of potential function. In each case we 
examine a first-order resonance, so the choice of a in the Lindstedt series is made when 
solving the order e equation. This also means the twist condition (17) is satisfied. 
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Figure 1. Logarithm of the jump in Dg^A^itj) against log(£') in the two 
dimensional example: V{x,u) - s sinilnkx) cosilnu), where k = (2,3). 
This is a "first order resonance" so w = fe. 256 modes were used for each 
Fourier frequency direction ^i, ^2- 




Figure 2. Plots of the logarithm of the jump in De,Ae(w) 
against log(e) for different potential functions. V(x, u) - 
esin(27rA;iXi)sin(27rA;2-«2)cos(27rM) on the left and V(x,u) - 
I sin(2;7rfcx)(cos(27rM) + sin(27rM)) on the right, k - (2, 1) in each 
case, and cj - k. 256 modes were used for each Fourier frequency 
direction ^i, ^2- 

Figure 1 is log-log plot of the jump in the e\ direction of the gradient of A£((y) versus e. 
The dotted line is the logarithms of the numerically computed points plotted against log(e). 
The soUd Une is the fit Jie) = \og{^e^^^), which was derived in Section 4.2. The same 
relation is found for the resonance at w = -k. 

Figure 2 has the same type of plot for different choices of potential function. In 2(a) the 
potential is 

y(x, u) = s sin(27r^ixi) sin(27r^2-^2) cos(2;7rM), 
and the choice of resonance was to = k - (2, 1). In this case, the fit is J(e) - log(^e'^^). 
The same relation is found at w = (+ki,+k2) and at w = (±^2. ±^i)- 
In 2(b) the potential is 

V{x, u) = - sin(27rA:x)(cos(27rw) + sin(2;rw)), 
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and the choice of to was a> - k - (2, 1). The same relation is found ato) - -k. 

Higher order resonances can be computed, for instance sXio -2k in each of the previous 
examples. However, the behavior of the jump in DAg{<j}) behaves like e to a power greater 
than one. To get reasonable accuracy one needs to take many more Fourier modes in the 
numerical approximation, which is impractical for us at the moment in the case of two 
spatial dimensions. 
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